arXiv:l506.04503v2 [astro-ph.GA] 17Jul2015 


Draft version July 20, 2015 

Preprint typeset using ETpX style emulateapj v. 6/22/04 


CALIBRATION OF THE OPTICAL MASS PROXY FOR CLUSTERS OF GALAXIES AND AN UPDATE OF THE 

WHL12 CLUSTER CATALOG 

Z. L. Wen, J. L. Han 

National Astronomical Observatories, Chinese Academy of Sciences, 20A Datun Road, Chaoyang District, Beijing 100012, China; zhonglue@nao.cas.cn. 

Draft version July 20, 2015 

ABSTRACT 

Accurately determining the mass of galaxy clusters is fundamental for many studies on cosmology and 
galaxy evolution. We collect and rescale the cluster masses of 1191 clusters of 0.05 < z < 0.75 estimated by 
X-ray or Sunyaev-Zeldovich measurements and use them to calibrate the optical mass proxy. The total r-band 
luminosity (in units of L*) of these clusters are obtained by using spectroscopic and photometric data of the 
Sloan Digital Sky Survey (SDSS). We find that the correlation between the cluster mass M 500 and total r-band 
luminosity L 500 significantly evolves with redshift. After correcting for the evolution, we define a new cluster 
richness Rl *,500 = £500 R(z) [ ' V) as the optical mass proxy. By using this newly defined richness and the 
recently released SDSS DR12 spectroscopic data, we update the WHL12 cluster catalog and identify 25,419 
new rich clusters at high redshift. In the SDSS spectroscopic survey region, about 89% of galaxy clusters have 
spectroscopic redshifts. The mass can be estimated with a scatter of 0.17 dex for the clusters in the updated 
catalog. 

Subject headings: galaxies: clusters: general — galaxies: distances and redshifts 


1. INTRODUCTION 

Clusters of galaxies are the most massive gravitationally 
bound systems in the universe. Statistics of cluster prop¬ 
erties provides very p owerful constrain t s on cosmology an d 
galaxy evolution (see I Allen et af] 1201 ll IWetzel et all 20121) . 
The cluster mass is a fundamental parameter for galaxy clus¬ 
ters and is related to observational f eatures in the X-ray, opti¬ 
cal, m illimeter and radio bands (e.g., Voit 2005; Bru netti et all 
120091) . Galaxies, intracluster hot gas and dark matter con¬ 
tribute about 5%, 15%, and 80% of the total cluster mass, re¬ 
spectively. Determining the cluster mass is the basis for many 
studies. For example, the mass function for a complete sam¬ 
ple of galaxy clusters can be used to constrain cosm ological 
parameters (e.g., Reinrich & Bohringer 2002 : Vikhlinin et ah 

l20^ht iWenet al .112010i). 

Statistics for weak gravitational lensing features is the most 
direct method to determine the total cluster mass including 
contributions from member galaxies, hot intracluster gas and 
dark matter. However, this method requires high-quality im¬ 
age data and also redshift information for faint background 
galaxies and thus only a few massive clusters have their 
masses so determi n ed (e.gJDahl e 2006; Barde au et alJl2007l : 
lOkabe et al.l 120101 : Ivon der Linden et akl l2014blk Based on 
the assumption that the hot intracluster gas is in hydro¬ 
static equilibrium, the cluster mass can be estimated by us¬ 
ing the X-ray surface brightness and temperature distributions 
(e.g.. iFinoguenov et alJ 12001b iReiprich & Bohringed 120021 : 

iRravtsov et all 120061 : IVikhlinin et al.l 120061) . A few thou¬ 
sands of galaxy clusters have their masses es t imated from 


the X-ray observables (IVikhlinin et al. 2009a| iMantz et~aH 
l2010t iPiffaretti et al.l 1201 It iTakev et al. 2013h . In the mil¬ 
limeter band, the thermal Sunyaev-Zeldovich (SZ) effect pro¬ 
vides robust ma ss estimates for more t han one thousand mas¬ 
sive clusters (lHas selfield et al.l 120131 : iReichardt et al.l 120131 : 
iPlanck C ollaborat ion et alJ 120141) . which is independent of 
cluster dynamical state and redshift (ICarlstrom et al.li2002l) . 

The velocity dispersion of member galaxies is certainly a 
good measure of cluster mass (iGirardi et al.lll998l) . which can 
be derived from optical spectroscopic data and have been ob¬ 


tained for a limited number of galaxy clusters. However, the 
spectroscopic observations are usually incomplete for clus¬ 
ter member galaxies, and clus ter substructures often induc e 
a bias on the estimated mass (lBirdlll995l : ISifon et aPl2013h . 
On the other hand, pho tometric data can be used to iden- 


tify galaxy clusters (e.g., Koester et a 

. 2007: Wen et al. 2009; 

Hao et al. 2010; Szabo et al. 2011; 

Wen et al. 

2012 Oguri 

2014; Rvkoff et al. 201-3). Cluster richness, defined as the 


total number of member galaxies or their total luminosity, 
can be used as an optical mass proxy (iPopesso et all 12005L 
120071 ). The challenge is to accurately determine the member¬ 
ship of cluster galaxies by eliminating the contamination by 
field galaxies. Based on multicolor survey data, the member 
galaxies can be discriminated by using galaxy colors (e.g. , 
iGladders & Ye3 l2000l lloester et alJ [20071 lHao et al.l 12010) 
or ph otometric redshifts (e.g.. IWen et alll2009h ISzabo et alJ 
1201 ill . The cluster richness derived from optical photometric 
data in general is poorly correlated with cluster mass. How¬ 
ever, if the correlation is improved, the cluster richness can 
be used to estimate cluster mass for a very large sample of 
clusters even up to high redshifts. __ 

The Sloan Digital Sky Survey (SDSS; [York et al .1 12001)1) of¬ 
fers an unprecedented photometric data in five broad bands 
(it, g, r, i, and z) covering 14,000 deg 2 with the excep¬ 
tional follow-up spectroscopic observations^ The photomet¬ 
ric data reach a limit of r = 22.2 (IStoughton et al.1 2002 ). 
with the star-galaxy separation reliable to a limit of r = 
21.5 (iLupton et aI. II 200 ft) . The spectroscopic survey observes 
galaxies with a Petrosian magn itude of r < 17.77 for the 
main galaxy sample (IStrauss et all 120021 and a model mag- 
nitude r < 19.9 for the Luminous Red Galaxy sample 
(lEisenstein et al.l 1200 ll) . These galaxies cover a footprint of 
10,400 deg 2 in the range of —11° < Decl. < 69°. The lat¬ 
est Data Release 12 (DR 12. lAlam et all 120151) includes the 
photometric data for about 200 million galaxies and the spec¬ 
troscopic data for about 2.32 million galaxies. 

The large samples of galaxy clusters or groups have 
been identified by using th e SPS S spectroscop i c data 
(e.g., iMerchan & Zandivarezl 120051: iBerlind et al.l 120061 : 
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Tempel et al. 

2014 

) and the photometric data (e.g., 

Koester et a 

2007 

; Wen et al. 2009; Hao et alJ 2010: 

Szabo et al. 2 

01 It IWen et all 2012 lOguril 20 ll Rvkoff et all 

2014[). Generally, the richness in these cluster catalogs 


is defined as the total number or the total luminosity of 
member galaxies, which can be regarded as a coarse mass 
proxy for clusters up to a redshift of z ~ 0.4. At higher 
redshifts, the discrimination of member galaxies are limited 
to only very luminous galaxies, and hence the richness is 
systematically underestimated, as discussed in IWen et ahl 
(120121) . We noticed that without a proper calibration, such 
a cluster ric hness could be redshift-dependent even at low 
redshift (e.g., IRozo et al.1120091 IWen et alKoIoh . 

Many recent samples of galaxy clusters have their 
masses determined by X-ray or SZ meas urements up 
to redshift of z ~ 0.75 (e.g., iPiffaretti et al.l 1201 It 

iPlanck Collaboration et alJl2014l) . These provide an oppor¬ 
tunity to calibrate the richness over a large redshift range so 
that the calibrated richness may be used as a mass proxy for 
a very large sample of galaxy clusters. In Section 2, we col¬ 
lect the cluster mass in the literature estimated by the X-ray 
or SZ observations, and for these clusters in the SDSS survey 
region, we calculate their total r-band luminosity and define 
a mass proxy that is independent of redshift. On the other 
hand, the SDSS Baryon Oscillation Spectroscopic Survey 
dDawson et al.ll2013h offers spectrosco pic redshif t s for most 
of the galaxy clusters in the catalog of lWen et al.l (1 20121 i.e., 
the WHL12 cluster catalog), so that in Section 3 we update 
the richness and redshift information for the WHL12 clusters, 
and find new complementary clusters at high redshift. 

Throughout this paper, we assume a ACDM cosmology, 
taking Ho = 100 h km s -1 Mpc -1 , with h = 0.7, fl m = 0.3 
and Ha = 0.7. 

2. CALIBRATING THE RICHNESS FOR AN OPTICAL MASS PROXY 

We collect the estimated masses of many large samples of 
galaxy clusters, and then scale them to each other by using 
common clusters. The scaled masses of the composite sample 
of galaxy clusters are used to calibrate the optical richness of 
galaxy clusters in a large range of redshifts so that any redshift 
dependence can be eliminated. 


2.1. The mass estimates for galaxy clusters and their 
rescaling between different samples 


Many samples of galaxy clusters have their mass e stim ated 
by us i ng the X-ray or SZ measurements (e.g.lKravtsov et all 
12005 17, hang et al 1120081: iPratt et al]12003 IZhao et. alll20m 

To make the largest composite sample with unified scal¬ 
ing on mass estimates, however, we collect only large sam¬ 
ples of clusters from the literature, with each having at least 
10 overlapping common clusters with the other samples so 
that the estimated masses of these common clusters can be 
compared and scaled. For this purpose, we take six large 


samples from t he liter ature: the X-ray cluster samples by 
IVikhlinin et al.l (l2009al hereafter V09 ), iMantz et all (120101 
hereafter M10), Piffaretti et~akl (1201 ll hereafter PI 1) and 
iTakev et akl (120131120141 hereafter T14), and the SZ clus¬ 
ter sample s from the Atacama Cos mology Telescope (ACT) 
survey by iHasselfielde^ak^ JgOjjj, hereafter H13) and the 
Planck survey by the lPlanck Collaboration et~akl d2015l here¬ 
after P15). We extract the cluster mass M 500 from these refer¬ 
ences, which is the mass within a radius of rsoo within which 
the mean density of a cluster is 500 times of the critical den¬ 


sity of the universe, p c = 3H 2 (z)/8nG, i.e. 

47T o 

M500 = Y r 500 x 500p c . (1) 

In the following we briefly introduce the six samples: 

V09: IVikhlinin et al.l (l2009al) estimated the total masses for 
85 clusters at z ~ 0.05 and z ~ 0.5 by using high- 
quality Chandra observations. The cluster mass M^oo was 
derived by using three X-ray proxies: average temperature 
Tx, ga s mass M gas , and Y x = T x x M gas . IVikhlinin et all 
(2009a|) calibrated their scaling relations with cluster mass 
(iKravtsov et al.1120061) by using a sample of well-observed re¬ 
laxed clusters at low redshifts. The mass estimates derived 
using the three proxies are well matched with very small scat¬ 
ters. We adopt the cluster mass M 500 they derived from M gas 
and the gas mass fraction, which depends weakly on cluster 
mass. 


M10: IMantz et~af.l ( 2010 ) used high-quality Chandra obser¬ 
vations to estimate the gas mass M gas for 130 massive clusters 
that have 0.1-2.4keV luminosity larger than 2.5 x 10 44 erg s _1 
in the redshift range of 0.05 < z < 0.88. Different from V09, 
they derived the total mass M 500 assuming a constant value 
of gas mass fraction. 


PI 1: IPiffaretti et aD(l201ll) compiled 1743 clusters in the red¬ 
shift range of z < 1.2 from t he ROSAT All Sky Su rvey- 
based d ata, including N ORAS dBohringe r et al.l I 2000I). RE ¬ 
FLE X dBohringer et alJ l200l . BCS (lEbeling et alJ l l 998h . 
SGP (ICruddace et alJl200l) . NEP dHenrv et al.li2006l). M ACS 
(lEbeling et all l200lT> and CIZA (lEbeling et alJ l2002h . and 
serendipito us searches in the ROSAT observations includ - 
ing 160SD (Mu llis et alJl2003l) . 400SD dBurenin et alJl200l 


SHARC (Rpmer et al 


and EMSS dGioia etlil 


2000i) WARPS ( Horner et al] 2008 1 
1990). All clusters have their X-ray 


luminosity, L x , 5 oo, measured within ?’ 5 oo- The cluster masses 
were estimated from L x , 500 by using the scaling relation be¬ 
tween mass and X-ray lumin osity obtained by Arnaud et al. 
(120101) and IPratt et all (120091) . with a typical uncertainty of 
15%. 


T14: ITakev et ahl (12013l) presented a sample of 530 X-ray 
clusters selected from the XMM-Newton serendipitous source 
catalog in the SDSS sky region, of which 310 clusters have 
spectroscopic redshifts and are in the redshift range of 0.03 < 
^ < 0.70. ITakev et alJ (120141) updated the photometric red¬ 
shifts of the clusters by using the SDSS DR10 spectroscopic 
data. Cluster masses are estimated by using the X-ray lumi¬ 
nosity L x 500 and the scaling relations given by IPratt et all 
(120091) . 


H13: lHasselfield et all (1201 3l) presented a catalog of 68 clus¬ 
ters over 0.15 < z < 0.8 detected via SZ effect at 148 GHz 
in the ACT survey. The authors provided cluster mass esti¬ 
mates for the 68 clusters near the celestial equator together 
with 23 clusters in the southern survey. Four different mass 
estimates were given in the paper based on the Compton pa- 
rameter calibrated to th e Universal Pressu re Profile model 
(lArnaud et al.ll2010l). thelBode et all (1201 2l) model, the Non- 
thermal mo del ( Trac et all201 ll) and the dynamical masses of 
ISifon et ah ( 2013 s). respectively. These four mass estimates 
are well correlated with different normalization. Here, we 
adopt the mass estimates calibrated to dynamical masses of 
ISif6n etal.l(f20T3h . 
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Fig. 1.— Comparison of the estimated masses for the overlapping common clusters in different samples from the literature. The solid line is the best fit, and 
the dotted line indicates the equivalence. All data in panels (a) — (e) are finally rescaled to V09 via the best fit(s), and summarized in panel (f), where the mass 
values in the X-axis are rescaled from the X-values of panels (a) - (e) and the m ass va lues in the Y-axis are rescaled from the Y-values of panels (a) - (e). Because 
the intrinsic uncertainty of mass-observable relation is larger than 5% jVikhlinin et al.l2009al) . we here set the minimum uncertainty of 5% for cluster mass even 
though it could be formally small in the literature. 


P15: The all-sky Planck catalog of clusters is the 

largest SZ-selected catalog to date, containing 1653 clus¬ 
ters with redshifts_up to z _~ 1 i n the second catalog 

JPlanck Collaboration et al .1120151) . Among these, 1203 clus¬ 
ters have counterparts in previous X-ray, optical, SZ, and in¬ 
frared data and 1094 clusters have redshifts. The masses 
of the 1094 clusters with redshift information were derived 
by using the Compton parameter dPlanck Collaboration et al.l 

12014 . 

The masses of galaxy clusters in these samples are com¬ 
piled to form a very large composite sample. However, we 
notice that cluster masses in these six samples were estimated 
by using different mass proxies with different scaling rela¬ 
tions, and even cluster masses derived from different X-ray 
proxies in the same sample are not consistent as well. Though 
both V09 and M10 derived cluster masses via the gas mass 
M gas , normalization of the fraction of gas mass and depen¬ 
dence on cluster mass are different in the two papers. Differ¬ 
ent cluster masses are obtained in HI3 if they are normalized 
to different models. Therefore, to get the uniformed scale for 
estimated masses of galaxy clusters in a composite sample, 
the mass estimates in different samples by different authors 
must be rescaled to each other for consistency. 

In general, the mass estimates derived by the proxies of 
M gas and Y\- are more reliable than those from the X- 
ray luminosity (iKravtsov et al.l 12006 *). For the SZ clus¬ 
ters, we find that the mass estimates depend on model 


for the ACT clusters dHasselfield et al.l 120131) ; the clus¬ 
ter masses from the Planck survey are probably sys¬ 
tematically underestimated ( iPlanck Collaboration et al 1120141: 
Ivon der Linden et aTll2014ai) . In this work, we rescale cluster 

masses to the standard obtained by IVikhlinin et al.l ( 2009a:) 
because their mass estimates are based on very high-quality 
X-ray data and the mass proxy shows only a small scatter af¬ 
ter careful calibration using the known masses of low redshift 
relaxed clusters. 

The easy approach to rescale the cluster masses to the stan¬ 
dard of lVikhlinin et al. (2009a ) is to compare the mass values 
of common clusters that overlap in these different samples. 
As shown in Figure Ufa), cluster masses derived by V09 and 
M10 using the same proxy M gas show a systematic offset. 
The best fit is found as 

logM5oo,vo9 = (0.99 ± 0.06) log Msoo.mio — (0.07 ±0.04). 

( 2 ) 

This equation can be used to convert the mass estimates of all 
galaxy clusters in the M10 sample to the V09 standard. Simi¬ 
larly, we compare the mass estimates for overlapping common 
galaxy clusters in V09 and P15, as shown in Figure [I] b), and 
obtain the conversion equation as 

logM5oo,vo9 = (1-06 ± 0.03) logM. 5 oo,pi 5 ± (0.02± 0.02). 

(3) 

In the cluster sample of PI 1, there are too few clusters that 
overlap with V09 sample. However, there are many overlap- 
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TABLE 1 

Rescaled mass and newly calculated richness for 1191 galaxy clusters 


ID 

R.A. 

Decl. 2 : 

M 500 

a -^500 

^500 

°> 5 oo 

Rl*, 500 

Reference 


(deg) 

(deg) 

(10 14 Mq) (10 14 Mq) (Mpc) (Mpc) 



(1) 

(2) 

(3) (4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

1 

0.49367 

12.06612 0.1995 

4.40 

0.42 

1.08 

0.04 

86.20 

P11.P15 

2 

0.79826 

-6.09169 0.2328 

7.20 

0.44 

1.26 

0.03 

105.02 

P11.P15 

3 

0.95698 

2.06647 0.0967 

2.57 

0.26 

0.94 

0.03 

60.13 

P11.P15 

4 

1.24353 

11.70091 0.0763 

2.00 

0.30 

0.87 

0.05 

23.18 

Pll 

5 

1.34984 

16.21926 0.1155 

3.70 

0.56 

1.05 

0.06 

67.62 

Pll 

6 

1.58453 

10.86429 0.1680 

4.46 

0.41 

1.10 

0.03 

68.22 

P11.P15 

7 

2.04332 

2.02009 0.3665 

6.44 

0.70 

1.16 

0.04 

68.70 

H13.P15 

8 

2.33609 

6.82309 0.2378 

5.80 

0.66 

1.17 

0.05 

112.20 

P15 

9 

2.51443 

17.77082 0.1714 

3.76 

0.59 

1.04 

0.06 

64.06 

P15 

10 2.57239 

6.64203 0.2648 

5.68 

0.72 

1.15 

0.05 

53.00 

P15 


Notes: Column (1): sequence number; Column (2) and (3): R.A. ('.12000J and Decl. (J2000) of cluster; Column (4): spectroscopic redshift of cluster; Column 
(5) and (6): rescaled cluster mass M 500 and its uncertainty o~m 500 ', Column (7) and (8): recalculated cluster radius r .500 and its uncertainty &r 50c from M 500 ; 
Column (9): newly defined cluster richness Rl*,500’ Column (10): reference for original mass estimate: V09 for IVikhlinin et al] 12009a 1 ). M10 for lMantz et all 
120101) . Pll for lPiffaretti et all i201 ID . H13 for lHasselfield et all<20131) . P15 for lPlanck Collaboration et aTl 120151) . T14 for lTakev et alj i20l3ll2014D . 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding the form 
and content.) 


ping common clusters with the P15 catalog. We can rescale 
the cluster masses of the Pll sample to the P15 scale first, 
and then convert to the V09 standard. Comparing the mass 
estimates of the common clusters overlapping in the Pll and 
P15 samples shown in FigureQJc) gives 

log Msoo.pis = (0.90 ± 0 . 02 ) log.M 500 .p 11 + (0.15 ± 0 . 01 ). 

(4) 

Combining the above two conversion equations, we get 

log M^5 oo,vo 9 = (0.95 ±0.03) log M 500 ,pii + (0.18 ±0.02). 

(5) 

The HI3 sample of the SZ clusters has even fewer overlaps 
with the samples of V09 and P15. However, there are some 
common clusters overlapping in the H13 and Pll samples. 
We compare their mass estimates as shown in Figure |TJd), 
and get 

logM 5 oo,ph = (0.98 ± 0.30) logMsoo.im — (0.19± 0.30). 

( 6 ) 

When it is combined with Equation (0, we get the conversion 
equation to the V09 standard as 

logM 50 o,vo9 = (0.93 ± 0.29) logM 50 o.Hi 3 ± (0-00 ± 0.29). 

(7) 

The same holds for the clusters in the T14 sample. The best 
fit between the T14 mass and the PI 1 mass as shown in Fig¬ 
ure He) is 

logM 5 oo,pn = (1.04 ± 0.08) log M 50 o,ti4 — (0.01 ±0.02), 

(8) 

and the conversion to the V09 standard is through 

logM5oo,vct9 = (0.99 ± 0.08) logM50o,Ti4 ± (0.17±0.03). 

(9) 

The masses of clusters thus scaled in Figure |T] a)-(e) have 
excellent consistency, as shown in Figure |TJ f), which shows 
a scatter of only 0.1 dex, i.e. <Ti og m 500 = 0.10. Most of 
the galaxy clusters from the V09, M10, Pll, H13, and P15 
samples have a rescaled mass of M 500 > 2 x 10 14 Mq, while 
clusters in the T14 sample often have a smaller mass of 0.3 x 
10 14 Mq < M 500 < 2 x 10 14 M e . 

The rescaled masses of the galaxy clusters from the six sam- 



FlG. 2.— The masses for 1191 clusters in the composite calibration sample 
are in the range of 0.3 - 20xl0 14 Mg, and the redshifts are distributed up 
to 0.75. 


pies will be used to calibrate the optical mass proxy, which 
will be derived from the SDSS data. Therefore we limit 
the composite sample to galaxy clusters in the SDSS sky re¬ 
gion, which must have a spectroscopic redshift either in the 
NASA/IPAC Extragalactic Database or from the SDSS spec¬ 
troscopic observations (lAlam et al.li2()15 ). We also limit the 
sample to a mass of M 500 > 0.3 x 10 14 Mq and in the red¬ 
shift range of 0.05 < z < 0.75. The lower redshift limit is 
set due to a photometry problem for very bright galaxies in 
the SDSS data, and the redshift upper limit is considered for 
the reliable detection of cluster galaxies in the SDSS data. 
We further visually inspect the color images 1 and exclude 
nine clusters that are seriously contaminated by very bright 
exposure-saturated stars. Finally we get 29, 64, 584, 50, 412, 
and 370 clusters from the V09, M10, Pll, H13, T14 and P15 
samples, respectively. After the overlapping common clus¬ 
ters are merged with their masses estimated from uncertainty- 
weighted average, the composite calibration sample consists 

1 http://skyserver.sdss. org/drl2/en/tools/chart/navi.aspx 
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Redshift, z 

Fig. 3.— Data scatter a z of photo -z varies with redshift, obtained from 
the comparison of spectroscopic redshifts and estimated photo -zs of a large 
SDSS galaxies. The dotted l ine r epresents the half thickness Az of the photo- 
z slice defined by Equation m for discrimination of member galaxies. 


of 1191 clusters, as listed in Table [I] The cluster radius r^oo 
are re-estimated from the rescaled mass M 500 by using Equa¬ 
tion 0. The distributions of their masses and redshifts are 
shown in Figure [2] 

2.2. Optical luminosity and richness of galaxy clusters 

We discriminate the cluster member galaxies and calculate 
the total optical luminosity of clusters by using the SDSS 
spectroscopic and photometric data following our previous 
procedures in tWen et all (120121 hereafter WHL12) but with 
a few improvements. The cluster richness is defined from the 
optical luminosity, which can act as the optical mass proxy 
for galaxy clusters. 


2.2.1. Discrimination of cluster member galaxies 

We only consider the bright member galaxies of M® ( 2 ) < 
—20.5. Here, M® is the evolution-corrected absolute magni¬ 
tude in the r band, M®(z) = M r (z) + Qz. Assuming the 
passive evolution of a galaxy popu lation formed at zy = 2 
(iLin et all l2006t ICrawford et al J 120091) and adopting a stellar 


population synthesis model ( Bruzual & Chariot! 120031 ) with 
the initial mass function of IChabrierl (120031) and solar metal- 
licity, we obtain the value of Q ~ 1.16 at redshift z < 0.8. 
The bright member galaxies of a cluster can be discriminated 
around a brightest cluster galaxy (BCG) in a small redshift 
range, i.e. a redshift slice. 

To discriminate the cluster galaxies, the spectroscopic red¬ 
shifts are taken if they are available in the SDSS DR12 
(iAlam et al.l:2QI 5): otherwise, photometric redshifts (photo- 
zs hereafter) will be used. For the objects classified as 
galaxies with only photo-zs, we remove those with deblend¬ 
ing problems and those marked as ‘SATURATED’ using the 
Hags 2 , and also discard those objects with a large photo-z er¬ 
ror zErr > 0.08(1 + z). These selection procedures remove 
most of the contamination from stars and galaxies with a bad 
photometry. 

For galaxies with spectroscopic redshifts, only those within 
a redshift slice of cAz/(l + z) = 2500kms _1 are considered 
as member galaxies. Here c is the speed of light, and z is the 
redshift of a cluster. 


2 (flags & 0x20) = 0 and (flags & 0x80000) = 0 and ((flags & 
0x400000000000) = 0 or psfmagenv <= 0.20) and ((flags & 0x40000) = 0) 


For most member galaxies, their photo-zs were estimated 
by optical magnitudes in five bands of the SDSS and have an 
uncertainty much larger than those of spectroscopic redshift 
(ICsabai et alJi2007l) . The data scatter of photo-z varies with 
redshift, as shown in Figure 0 and is smaller than 0.035 at 
a redshift of z < 0.5 but increases to a z ~ 0.075 at a red¬ 
shift of z > 0.7. Only galaxies within the redshift slice of 
z ± Az are con sidered as m ember gal axies of a c luster at a 
redshift of z. In I Wen et alJ (l2009h and IWen et al.1 (1201 2l) the 
half thickness of the slice is taken as Az = 0.04 (1 + z), and 
the majority of the member galaxies can be well discriminated 
for clusters with a redshift of z < 0.42. However, due to the 
large uncertainty of photo-z, such a slice leads to incomplete 
recognizance of member galaxies for clusters at higher red¬ 
shifts. Here we adopt the photo-z slice in the form of 


f 0.04 (1 + z) for z< 0.45 
\ 0.248 z- 0.0536 for z> 0.45 ’ " 


which is about 1.75-2.00 a z and by which most cluster mem¬ 
ber galaxies can be identified up to a high redshift. Such a 
larger photo-z slice leads to a slightly higher completeness for 
identification of member galaxies but introduces more con¬ 
tamination. 

Compared to IWen et aD (2012), we improve the discrimi¬ 
nation of member galaxies in three aspects: (a) spectroscopic 
redshifts are considered for member galaxies; (b) the evolu¬ 
tion of abso lut e magnitu des is improved from Q = 1.62 ± 
0.30 (Blanton et al.ll2003 0 to the current value of Q = 1.16; 
and (c) the thickness of the photo-z slice at high redshifts is 
enlarged according to the photo-z scatter. 


2.2.2. Total r-band luminosity of cluster member galaxies 
and background estimation 

In general for each known cluster, the total r-band lumi¬ 
nosity can be calculated by summing the luminosity of rec¬ 
ognized member galaxies. The total luminosity is measured 
in units of L* which is the evolved characteristic luminosity 
of galaxies in the r-band defined by L*(z) = L*(0)10°’ 4< 5 Z , 
here again Q = 1.16. 

However, the contamination by background field galaxies 
has to be subtracted. To estimate the local bac kgroun d, we 
follow the method analogous to iPopesso et al.1 ( 2004 ). For 
each cluster, we take the BCG as the center of a galaxy clus¬ 
ter, and then divide the annulus between 2-4 Mpc from the 
BCG into 48 sectors with equal areas. Within the same mag¬ 
nitude limit and the same photo-z slice, the field galaxies in 
these sectors fainter than the second brightest cluster member 
galaxies are used to calculate the background. We get the total 
r-band luminosity of bright galaxies in each sector, and esti¬ 
mate the mean and deviation a. The sectors with a luminosity 
obviously deviating from the mean by more than 3cr are dis¬ 
carded, and the mean is recalculated for the local background 
which is scaled and subtracted from the total luminosity. 


2.2.3. Scaling relations between the optical data and mass 

The goal here is to establish the scaling relation between the 
cluster mass and the luminosity or richness derived from op¬ 
tical data. As shown bv lWen et al.1 (. 2012 "). if there are merely 
optical survey data for a cluster, a few steps are needed to de¬ 
rive the richness from optical data: (1) discriminate cluster 
member galaxies and get their redshifts; (2) get the optical 
luminosity Lim P c from member galaxies within 1 Mpc, and 
then derive the cluster radius rsoo; and (3) get the optical lu- 
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log E(z) log E(z) log E(z) 

Fig. 4.— Upper panels are the scaling relations between cluster parameters derived from the optical data (the X-axis) and cluster radius or mass (the Y-axis) 
for the 1191 clusters from the composite calibration sample. The data scatters are plotted in the lower panels against log E(z) to check their possible evolution. 


minosity ±500 within 7 - 500 , and then the richness defined as 

Rl *,500 = L 500 /L*. 

First, to estimate rsoo from optical data, we have to estab¬ 
lish the relation between r $00 and an optical observable that is 
the total r-band luminosity L \ m pc (in units of ±*, with back¬ 
ground subtracted) of member galaxies within 1 Mpc from the 
cluster center. As shown in the left panels of Figure[4] we find 
that the values of rsoo (in Mpc) in Tableware well correlated 
with -LiMpc- The uncertainty of ±!M pc comes from the in¬ 
completeness of member galaxies and residual contamination 
of background galaxies, which is about 10% for rich clusters 
dWen et al.ll2009l ). By taking a typical uncertainty of 15% on 
TiMpc f°r all clusters, we find the best fit to the data as 

logrsoo = (0.45 ±0.01) log ±im pc - (0.76 ±0.01). (11) 

The slope obtained here is consistent with that in IWen et akl 
(120121) . 

Some authors have shown that the scaling relation between 
Tf .500 (°r r. 5 oo) and an optical observable (e.g., richnes s) may 
evolve with redshift (|Lin et alJ 120041 lAndreon & Co ngdon 
l2014t IWiesner et al.ll2015h , for example, due to higher con¬ 
centration of the light-to-mass profiles at lower redshifts 
iGu et aljl2013l) . On the other hand, bright member galaxies 
of M e < —20.5 are incomplete in the SDSS photometric data 
at z > 0.45, so that the total luminosity could be underesti¬ 
mated. The 1191 clusters in the composite calibration sample 
are in a very large redshift range, allowing us to clarify any 
redshift evolution of the scaling relations. For the scaling re¬ 
lation rsoo-±iM P c> we calculate the deviation of r $00 from the 
relation: Alogr 500 = logr 500 -(0.45 log ±im P c -0.76) and 
find that A log r 500 is only slightly evolving with the redshift 
in the form of (see the lower left panel in Figure [4} 

A logrsoo = (-0.20±0.06)log±;(z)±(0.01±0.01), (12) 

where E(z) = -^/Ha ± fl m (l + z) 3 . Including this small 
correction, the radius rsoo is related to L \ m pc and z by 

logrsoo = (0.45 ± 0.01) log± 1Mp c - (0.75 ± 0.01) 

— (0.20 ± 0.06) logE(z). (13) 


We calculate the total r-band luminosity ±500 (also in 
units of L*, with background subtracted) of member galax¬ 
ies within the radius of rsoo for the 1191 clusters by using the 
SDSS data. As shown in the upper middle panel of Figure [4] 
there exists excellent correlation between ±500 and the cluster 
mass M 500 (in units of 10 14 M 0 ), and the best fit is 

logMsoo = (1.08 ± 0.02) logL 500 - (1-28 ± 0.02). (14) 

Again, to check for the possible evolution with redshift, 
the deviations of M 500 from the above fitted relation, 
AlogM .500 = log M 500 - (1.08 log L 500 - 1-28), are plot¬ 
ted against E(z) in the middle lower panel of FigureQ] which 
shows a correlation in the form of 

A log M 500 = (1.51 ±0.14) log±(z) — (0.09 ±0.01), (15) 

Though the total r-band luminosity ±500 in principle may 
slightly be underestimated for clusters of z > 0.45 

(i.e., log ±(0.45) = 0.104), the A logMsoo has to be so cor¬ 
rected. We then can define the richness of the galaxy clus¬ 
ter from the total r-band luminosity ± 500 , which can act as a 
well-calibrated mass proxy in the optical band, as 

Rl*,500 = L 500 E(z) 1A0 . (16) 

Such a richness is independent of redshift and is valid in the 
redshift range of 0.05 < z < 0.75, as shown in the right 
panels of Figure[4] The thus-calculated richness Rl*, 500 for 
galaxy clusters in the composite calibration sample is listed in 
TableQ] The scaling relation between the cluster mass and the 
richness is 

logMsoo = (1.08±0.02)logi?A* j5 oo — (1.37± 0 . 02 ), (17) 

which is consistent with but more accurate than the mass- 
richn ess relation obtained in our previ ous work (IWen et alJ 
120121) and the M 500 -A relation given bv lRvkoff et alJ (120121) . 
The scatter of mass values is about 0.17 dex, i.e. <j\ O gM 500 = 
0.17. Taking out the inherent scatter (J\ og M 500 = 0.10 of the 
rescaled masses of the composite calibration sample, the mass 
uncertainty estimated by Rl*, 500 should be <Ji og M 500 = 0.14. 

Because clusters in the composite calibration sample are 
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log E(z) 

Fig. 5.— Correlation between M 500 and the directly optically estimated 
^L*, 5 QQ for the 1191 clusters in Table |T| Different from the right panel of 
Figure \4\ the values of 7-500 and Rl * 500 are both estimated from merely 
optical data. 


distributed in a wide range of mass and redshift, the scaling 
relations we derived above should be valid for clusters in the 
SDSS in the redshift ranges of 0.05 < z < 0.75 with a mass 

of M 500 > 0.3 x 10 14 Mq. 

2.3. Goodness of optical mass proxies 

How good is the richness Rl *,500 estimated from just the 
optical data (i.e., 7-500 is unknown previously) acting as a mass 
proxy? To get a quantitative estimate, we estimated the L\ m pc 
of 1191 clusters in the composite calibration sample, and then 
get 7-500 from the scaling relation in Equation ([13} so that 
L 500 and hence Rl *,500 can be estimated from the member 
galaxies within 7 - 500 - Obviously the uncertainty of estimated 
7-500 is propagated to L 500 and then Rl*, 500 - The plot for 
M 500 and Rl *,500 in Figure [5] shows a mass scatter in 0.20 
dex, i.e. a\ og M 50 o = 0.20. Taking out the inherent scatter of 
< 7 tog M 500 = 0.10 for the masses estimated from the X-ray and 
SZ observations, the mass scaled by Rl *,500 should have a 
scatter of cr logM500 = 0.17. 

Can the richness in other cluster catalogs be used as a 
mass proxy? If so, how good is it? To find out the an¬ 
swer to these questions, we investigate the optical richness 
based on photometric data in the cluster catalogs of the 
maxBCG dKoester et al.ll2007h . the Gaussian Mixture Bright¬ 
est Cluster Galaxy (GMBCG, lHao et all 1201 Oh . the adap¬ 
tive matched fil ter (AMF, ISzabo et all 201 lh . the WHL12 
(IWen et al.l2012h . the red-se quence Matched-filt er Probabilis¬ 
tic Percolation (redMapper, Rvkoff et al.ll2014h . the Cluster 
finding Algorithm based on Multi-band Identification of Red- 
sequence galaxies (CAMIRA, lOguril 120141). an d the rich¬ 
ness based on spectroscopic data by lYang et all ( 2007 ) and 


iTempel et al.l (120141) . 

The richness is defined as the total number of galaxies in the 
catalogs of maxBCG, GMBCG, and redMapper, or defined as 
the total luminosity of member galaxies in the catalogs of the 
AMF and WHL12. lYang et all (120071) adopted the charac¬ 
teristic luminosities or the characteristic stellar masses as the 
halo mass proxy of galaxy group/clusters, and here we take 
the halo mass for comparison. iTempel et al.l (12014h provided 
the dynamical mass estimates based on radial velocity disper¬ 
sion of member galaxies, which we will take for comparison. 

We cross-match the galaxy clusters in Table 1 with clusters 
in the above cluster catalogs, and get the matched clusters 
within a projected distance of 7-500 and a redshift difference 
of 0.05. If more than one cluster in a catalog matches with 
a cluster in Table 1, we take the one with the minimum pro¬ 
jected distance. Note that for different cluster catalogs, we 
take the redshift ranges up to their redshift distribution peaks. 
As shown in Table [2] we get different numbers of matched 
clusters for different catalogs. We estimate the cluster mass 
by using each richness and get the scatter of the mass. For 
comparison, the mass scatter for the same sample of matched 
clusters for each catalog is also calculated by using our newly 
defined richness. We find that the richness defined in this pa¬ 
per and those in redMapper and CAMIRA are equally good 
to act as optical mass proxy. The improvement in this pa¬ 
per is that the newly defined richness has the redshift evolu¬ 
tion corrected and is good for a wide range of redshifts up to 
z - 0.75. 

3. UPDATE OF THE WHL12 CLUSTER CATALOG 

We identified 132,684 galaxy clusters by using the photo- 
z of galaxies from SDSS DR8 (IWen et al.l 1201 21) . Clusters 
were identified if they have a richness > 12 and the number 
of member galaxies A^oo > 8 within a radius of 7 - 200 - Here, 
7-200 is the radius within which the mean density of a cluster 
is 200 times of the critical density of the universe. The cluster 
richness Rl* was defined to be the r-band total luminosity 
(in units of L*) of member galaxies with M® < —20.5 within 
f200- 

Two reasons lead us to update the WHL12 cluster catalog. 
First, the latest SDSS DR12 (lAlam et al .1120 1 5l) released spec¬ 
troscopic redshifts for about 2.30 million galaxies, of which 
about 1.35 million are luminous red galaxies (LRGs). In 
general LRGs are massive galaxies, and many of them are 
the BCGs or the bright member galaxies of clusters in the 
WHL12. We can update these spectroscopic redshifts for 
most of the WHL12 clusters. Second, the newly defined rich¬ 
ness can act as an excellent mass proxy, and should be calcu¬ 
lated for all clusters in the WHL12 catalog. In addition, by 
using the SDSS DR12 data, we can find rich galaxy clusters 
at high redshifts. 

3.1. Update on spectroscopic redshifts and the BCGs 

In general, we take the spectroscopic redshift of the BCG 
t's.BCG as the redshift of a cluster. SDSS DR12 h as sp ectro¬ 
sco pic redshifts for about 2.30 million galaxies (lAlam et aD 
120151) . Simple cross-matching the BCGs of the WHL12 clus¬ 
ters with the galaxies with spectroscopic data in the SDSS 
DR12 gives the redshifts for 84,760 (64%) clusters. 

We then try to get spectroscopic redshifts of bright member 
galaxies. For clusters with available z s ,bcg> we consider only 
galaxies within 7-500 = 2/3 7-200 from the BCG and within 
2500 km s _1 from z s bcg in the rest frame. Here 7-200 can 
be found in Table 1 of tWen et al.l 120121) . For clusters without 
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TABLE 2 

Comparison of data scatter for masses estimated from different richnesses in different cluster catalogs. 

Cluster catalog Redshift range Number of <7i og m 500 CT iog m 500 

for matching matched clusters scaled from ref. from this work 


maxBCG 0.10 < 

z < 0.30 

325 

0.29 

0.20 

GMBCG 0.10 < 

z < 0.37 

502 

0.27 

0.19 

AMF 0.05 < 

z < 0.38 

402 

0.26 

0.20 

WHL 12 0.05 < 

z < 0.42 

799 

0.24 

0.20 

redMapper 0.08 < 

z < 0.55 

610 

0.17 

0.18 

CAMIRA 0.1 < z < 0.6 

647 

0.18 

0.19 

Yang et al. 

z < 0.2 

195 

0.56 

0.19 

Tempel et al. 

z < 0.2 

149 

0.23 

0.19 

This paper 0.05 < 

z < 0.75 

1191 

- 

0.20 



Fig. 6.— Distribution of the newly defined richness Rl* 500 for the 
WHL12 clusters. 


~s,bcg> we consider the galaxies within 7*500 =2/3 7*200 from 
the BCGs and within a redshift difference of 0.03(1 + z) from 
z p , the photo-z of clusters in IWen et all (120121) . This small 
redshift difference is about 2.0-2.5 times the uncertainty of 
the cluster photo-z. The mean spectroscopic redshifts of the 
member galaxies (including the BCG) is regarded to be the 
spectroscopic redshift of a cluster. Through this process, we 
get spectroscopic redshifts for 95,684 clusters, which is about 
72% in the WHL catalog. In the SDSS spectroscopic survey 
region, 85% of the WHL clusters have spectroscopic redshifts 
obtained. Clusters located at low Galactic latitudes have no 
spectroscopic redshifts yet. 

The position and the BCG have been updated for some clus¬ 
ters as well. In lWen et al.ld2012l) , all objects with the deblend¬ 
ing pro blem s or marked as “SATURATED” (see the flags in 
Section 12.2b were removed from the photometric data in the 
procedure for identifying galaxy clusters. The true BCGs that 
have such a flag or catastrophic photo-z have now been in¬ 
cluded by using the spectroscopic redshifts. If a galaxy is 
brighter than the BCG of an original WHL12 cluster within 
7*200 and has a redshift difference of 0.03(1 + z) and is sur¬ 
rounded by other member galaxies in the SDSS color image, 
we take it to be the new BCG of the cluster. We hence update 
the BCGs for 3128 clusters. 

3.2. Update on cluster richness 



Redshift, z 

Fig. 7.— Redshift distribution of the newly identified clusters, compared 
with the WHL 12 clusters of Rl* 500 > 20 in the SDSS spectroscopic survey 
region. 

By using the scaling relations we obtained in Section [2.2.31 
we update the richness of galaxy clusters in the WHL12 cata¬ 
log via the following steps: 

ilMpc —> T500 —> i.500 —I Rl *,500 ■ 

That is to say, we have to estimate r 500 from LiMpc. and then 
get L 500 and Rl*, 500 - 

The updated parameters are listed in Table [3j where the 
cluster names, coordinates of the BCG, cluster redshift, and 
the r-band magnitude, as well as the newly estimated radius 
7*500 and newly defined richness -Rl *,500 are listed. The num¬ 
ber of member galaxies JV 500 within r 50 o and the number of 
these galaxies with spectroscopic redshifts JV 5 qq sp are given 
for reference. The distribution of Rl *,500 is shown in Fig¬ 
ure [6j it has a peak around Rl *,500 = 20. We noticed that 
5183 clusters in the WHL12 catalog has a low new richness 
of Rl *,500 < 8. 

3.3. Complementary galaxy clusters at high redshifts 

Identification of galaxy clusters at high redshifts was in¬ 
complete due to the large uncertainty of the photometric red¬ 
shift, so that the the number of clusters in the WHL 12 ca talog 
drops down at redshifts of z > 0.42 dWen et alJl20Ll see 
also Figure [7]l. Here we identify complementary rich clus¬ 
ters around bright galaxies which have spectroscopic redshifts 
measured in the SDSS DR12, with following procedures: 
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TABLE 3 

The WHL12 cluster catalog with updated parameters 


Cluster name 

R.A. 

Decl. 

Zp 


iBCG 

r200 

Rl* 

(V200 

rsoo 

Rl*, 500 

(Woo.sp 

(V500 


(deg) 

(deg) 




(Mpc) 



(Mpc) 




(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(?) 

(8) 

(9) 

(10) 

(ID 

(12) 

(13) 

WHL J000000.6+321233 

0.00236 

32.20925 

0.1274 

-1.0000 

14.92 

1.72 

70.63 

24 

1.01 

68.28 

0 

19 

WHL J000002.3+051718 

0.00957 

5.28827 

0.1696 

0.1694 

16.20 

0.94 

17.48 

9 

0.72 

24.46 

1 

9 

WHL J000003.1—033245* 

0.01276 

-3.54578 

0.5968 

0.5973 

20.14 

1.00 

19.19 

11 

0.83 

49.33 

2 

13 

WHL J000003.3+311354 

0.01377 

31.23175 

0.5428 

0.5293 

20.17 

0.87 

14.27 

8 

0.63 

25.43 

2 

9 

WHL J000003.5+314708 

0.01475 

31.78564 

0.0932 

0.0916 

15.18 

0.94 

16.97 

9 

0.69 

21.14 

1 

11 

WHL J000004.7+022826 

0.01945 

2.47386 

0.4179 

-1.0000 

19.32 

0.95 

13.71 

10 

0.54 

6.60 

0 

3 

WHL J000005.5+354610 

0.02303 

35.76957 

0.4762 

-1.0000 

19.59 

0.91 

15.58 

9 

0.63 

19.13 

0 

6 

WHL J000006.0+152548 

0.02482 

15.42990 

0.1656 

0.1731 

16.60 

1.13 

23.53 

19 

0.79 

19.53 

2 

12 

WHL J000006.3+221220 

0.02643 

22.20558 

0.3985 

-1.0000 

19.36 

0.84 

12.73 

11 

0.52 

11.92 

0 

7 

WHL J000006.6+100648 

0.02755 

10.11333 

0.3676 

0.3747 

19.07 

0.93 

16.73 

13 

0.70 

23.52 

2 

12 


Note. Column (1): Cluster name, the name with **’ means the BCG is updated; Column (2) and (3): R.A. (J2000) and Decl. (J2000) of cluster; Column (4): 
photometric redshift of cluster; Column (5): spectroscopic redshift of cluster, —1.0000 means not available; Column (6): r-band magnitude of the BCG; Column 
(7): original r 2 oo of cluster (Mpc); Column (8): original richness; Column (9): original number of member galaxy within 7 * 200 ; Column (10): cluster radius 7*500 
derived from optical luminosity; Column (11): newly defined cluster richness, Rl* ,5001 Column (12): number of member galaxies with spectroscopic redshifts 
within 7*500; Column (13): number of all member galaxies within 7 * 500 - Column (4) and (7)-(9) are from lWen et alJ J2012D . Column (10)—(13) are derived in this 
paper. 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding the form and content.) 


TABLE4 

Newly identified rich clusters at high redshifts 


Cluster name 

R.A. 

Decl. 

z 

r*BCG 

r 500 

Hi *,500 

(V500,sp 

(V500 


(deg) 

(deg) 



(Mpc) 




(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

WH J000004.2+021941 

0.01743 

2.32800 

0.6443 

20.15 

0.69 

27.88 

1 

6 

WH J000004.3—091120 

0.01809 

-9.18902 

0.6028 

19.95 

0.70 

34.93 

1 

7 

WH J000005.9+100144 

0.02452 

10.02896 

0.6712 

20.46 

0.65 

27.84 

1 

6 

WH J000006.8+195244 

0.02824 

19.87897 

0.4771 

19.88 

0.62 

22.75 

2 

7 

WH J000007.7+234443 

0.03208 

23.74541 

0.5285 

19.84 

0.65 

27.25 

2 

7 

WH J000008.2+263228 

0.03423 

26.54098 

0.5462 

20.02 

0.64 

26.35 

1 

8 

WH J000008.6+053257 

0.03568 

5.54922 

0.5742 

20.05 

0.73 

34.66 

3 

10 

WH J000009.9—072520 

0.04132 

-7.42227 

0.5847 

20.78 

0.62 

24.19 

1 

8 

WH J000010.4+041038 

0.04318 

4.17714 

0.6954 

21.29 

0.58 

24.96 

1 

6 

WH J000013.5—022138 

0.05613 

-2.36069 

0.4959 

19.94 

0.61 

22.55 

2 

7 


Note. Column (1): cluster name; Column (2) and (3): R.A. (J2000) and Decl. (J2000) of cluster; Column (4): spectroscopic redshift of cluster; Column (5): 
r-band magnitude of the BCG; Column (6): cluster radius 7*500 derived from optical luminosity L\ Mpc5 Column (7): cluster richness, Rl*, 500 , Column (8): 
number of member galaxies with spectroscopic redshifts within 7 * 500 ; Column (9): number of all member galaxies within 7 * 500 - 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding the form and content.) 


1. Get the spectroscopic redshift of a cluster candidate. 
Each bright galaxy of M e < —20.5 with a spectroscopic red¬ 
shift is assumed to be a member galaxy of a cluster candidate. 
All galaxies within 0.5 Mpc from a given galaxy and with a 
velocity difference in the rest frame less than 2500 km s -1 
are taken to be the member galaxies of the cluster candidate. 
Then the redshift of the cluster is defined as the median value 
of the spectroscopic redshifts of these galaxies. 

2. Determine the BCG. Generally, the true BCG, because 
it is the brightest, should have been observed in the spectro¬ 
scopic survey already. However, the BCGs in some clusters 
may be not observed spectroscopically because of their blue 
colors (ICrawford et alJ l!999: Pip ing et al.l 1201 ll) or the fiber 
collision limitation in the SDSS spectroscopic survey. The 
BCG candidates are the galaxies within 0.5 Mpc and with 
a velocity difference less than 2500 km s -1 from the target 
galaxy or with a photo-:; difference of 0.03(1 + z). The small 
slice of photo-;: is used here for the BCG selection to reduce 
contamination from field galaxies. The BCG of a cluster can¬ 


didate is recognized as the brightest galaxy of these galaxies. 

3. Calculate the richness and judge a rich cluster. After we 
have the redshift of a cluster candidate and also the recognized 
BCG, we then find the bright member galaxies in the r edshift 
slice defined by Equation [lO] as described in Section 12.2.31 
and then calculate CiMpc to estimate r 500 , and then get L 500 
and finally Rl *, 500- We only identify rich clusters with a rich¬ 
ness of Rl *,500 > 20 (beyond the peak in Figure [ 6 } and the 
number of bright galaxies within J’soo- W500 > 6. 

4. Clean the entries. It is possible that one rich cluster 
can be recognized two or more times in the above procedures 
when targeting different bright galaxies, and hence some en¬ 
tries should be merged to one cluster if their redshift differ¬ 
ence is less than 0.1 and a projected distance is less than 

1-5 rsoo- 

We find 79,498 clusters through the above procedures, of 
which 54,079 can match the clusters in the WHL12 catalog 
with a redshift difference of 0.1 and a projected distance of 
1.5 '/"sou- The remaining 25,419 are newly identified clusters. 
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Fig. 8.— Sky coverage of rich clusters of Rl *,500 > 20 (upper panel) 
and their radial distribution in a slice 0° < Decl. < 5° (lower panel). The 
red points are clusters with spectroscopic redshifts. The dot size in the lower 
panel is scaled by ^/Rl*, 500 - 


as listed in Table Q] Most of them have a high redshift of 
2 > 0.4 as shown in Figure[7] 

If we combine the cluster sample in the updated WHL12 
catalog and the complementary catalog of newly identified 
25,419 clusters, we get 158,103 clusters, of which 96,257 
have a richness of Rl *,500 > 20 . Figure [ 8 ] shows the sky 
coverage and radial distribution in a sky slice for clusters in 
the combined catalog. In the SDSS spectroscopic survey re¬ 
gion, 89% of clusters have spectroscopic redshifts. 

4. CONCLUSIONS 

In this paper, we calibrated the optical mass proxy of galaxy 
clusters and applied it to update the WHL12 cluster catalog. 

We collected clusters with known masses (M 500 ) derived 
from X-ray or SZ measurements from the literature, and 
rescaled them for consistency. The composite mass calibra¬ 
tion sample consists of 1191 clusters in the redshift range 
0.05 < z < 0.75 with mass M 500 > 0.3 x 10 14 Mq, and 


all of them are limited to the SDSS sky coverage to obtain the 
optical spectroscopic and photometric data. 

We established the scaling relations between the total r- 
band optical luminosity of member galaxies L iMpc within 
1 Mpc and the cluster radius (rsoo). After correcting for 
weak redshift dependence, the scaling relation can be used 
to estimate the cluster radius from optical data. The cluster 
mass M 500 is found to be well-correlated with the total r- 
band luminosity L 500 , which evolves significantly with red¬ 
shift. Considering the redshift evolution, we defined the new 
cluster richness as Rl*, 5 00 = £500 E(z) 1A0 , which can be 
used as a mass proxy independent of redshift and is valid in 
the redshift range of 0.05 < z < 0.75. The scatter of mass 
estimates is about 0.17 dex. 

We updated the spectroscopic redshifts for the WHL12 
clusters using the SDSS DR12 spectroscopic data, and then 
obtain their newly defined richness. We identify 25,419 com¬ 
plementary clusters at high redshifts around bright galaxies 
with a spectroscopic redshift in the SDSS DR12. The finial 
combined catalog contains 158,103 clusters in total, and about 
89% of them have spectroscopic redshifts in the SDSS spec¬ 
troscopic survey region. 
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